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In this paper, we consider smooth shot noise processes and their 
expected number of level crossings. When the kernel response func- 
tion is sufficiently smooth, the mean number of crossings function 
is obtained through an integral formula. Moreover, as the intensity 
increases, or equivalently, as the number of shots becomes larger, a 
normal convergence to the classical Rice's formula for Gaussian pro- 
cesses is obtained. The Gaussian kernel function, that corresponds to 
many applications in physics, is studied in detail and two different 
regimes are exhibited. 

1. Introduction. In this paper, we will consider a shot noise process 
which is a real-valued random process given by 

(l) X(t) = J2Pi9(t-Ti), teR, 

i 

where g is a given (deterministic) measurable function (it will be called 
the kernel function of the shot noise process), the {t{\ are the points of a 
Poisson point process on the line of intensity Av(ds), where A > and v is 
a positive u-finite measure on K. and the are independent copies of a 
random variable (3 (called the impulse), independent of {t{}. 

Shot noise processes are related to many problems in physics as they result 
from the superposition of "shot effects" which occur at random. Fundamen- 
tal results were obtained by Rice [23]. Daley [10] gave sufficient conditions 
on the kernel function to ensure the convergence of the formal series in a 
preliminary work. General results, including sample paths properties, were 
given by Rosihski [24] in a more general setting. In most of the literature the 
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measure v is the Lebesgue measure on R such that the shot noise process is a 
stationary one. In order to derive more precise sample paths properties and 
especially crossings rates, mainly two properties have been extensively exhib- 
ited and used. The first one is the Markov property, which is valid, choosing 
a noncontinuous positive causal kernel function, that is, for negative time. 
This is the case, in particular, of the exponential kernel g(t) = e~ l tt>o for 
which explicit distributions and crossings rates can be obtained [21]. A sim- 
ple formula for the expected numbers of level crossings is valid for more 
general kernels of this type but resulting shot noise processes are nondif- 
ferentiable [4, 16]. The infinitely divisible property is the second main tool. 
Actually, this allows us to establish convergence to a Gaussian process as 
the intensity increases [15, 22]. Sample paths properties of Gaussian pro- 
cesses have been extensively studied and fine results are known concerning 
the level crossings of smooth Gaussian processes (see [2, 9], e.g.). 

The goal of the paper is to study the crossings of a shot noise process in 
the general case when the kernel function g is smooth. In this setting we lose 
Markov's property but the shot noise process inherits smoothness properties. 
Integral formulas for the number of level crossings of smooth processes was 
generalized to the non-Gaussian case by [18] but it uses assumptions that 
rely on properties of some densities, which may not be valid for shot noise 
processes. We derive integral formulas for the mean number of crossings 
function and pay a special interest in the continuity of this function with 
respect to the level. Exploiting further on normal convergence, we exhibit 
a Gaussian regime for the mean number of crossings function when the 
intensity goes to infinity. A particular example, which is studied in detail, 
concerns the shot noise process where /3 = 1 almost surely and g is a Gaussian 
kernel of width a, 



Such a model has many applications because it is solution of the heat equa- 
tion (we consider a as a variable) , and it thus models a diffusion from random 
sources (the points of the Poisson point process). 

The paper is organized as follows. In Section 2, we consider crossings 
for general smooth processes. We give an explicit formula for the Fourier 
transform of the mean number of crossings function of a process X in terms 
of the characteristic function of (X(t), X'(t)). One of the difficulties is then 
to obtain results for the mean number of crossings of a given level a and 
not only for almost every a. Thus we focus on the continuity property of 
the mean number of crossings function. Section 3 is devoted to crossings for 
a smooth shot noise process X defined by (1). In order to get the continuity 
of the mean number of crossings function, we study the question of the 
existence and the boundedness of a probability density for X(t). In Section 4, 



g{t)=9a{t) 
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we show how, and in which sense, the mean number of crossings function 
converges to the one of a Gaussian process when the intensity A goes to 
infinity. We give rates of this convergence. Finally, in Section 5, we study in 
detail the case of a Gaussian kernel of width a. We are mainly interested 
in the mean number of local extrema of this process, as a function of a. 
Thanks to the heat equation, and also to scaling properties between a and 
A, we prove that the mean number of local extrema is a decreasing function 
of a, and give its asymptotics as a is small or large. 

2. Crossings of smooth processes. The goal of this section is to investi- 
gate crossings of general smooth processes in order to get results for smooth 
shot noise processes. This is a very different situation from the one studied 
in [[4], [16, 21]] where shot noise processes are nondiffcrcnt iablc. However, 
crossings of smooth processes have been extensively studied, especially in the 
Gaussian processes realm (see [2], e.g.) which are second order processes. 
Therefore, in the whole section, we will consider second order processes 
which are both almost surely and mean square continuously differentiable 
(see [1], Section 2.2, e.g.). This implies, in particular, that the derivatives 
are also second order processes. Moreover, most of known results on cross- 
ings are based on assumptions on density probabilities, which are not well 
adapted for shot noise processes. In this section, we revisit these results with 
a more adapted point of view based on characteristic functions. 

When X is an almost surely continuously differentiable process on ]R, we 
can consider its multiplicity function on an interval [a, b] defined by 

(2) V«£l N x (a,[a,b])=#{t<E[a,b];X{t) = a}. 

This defines a positive random process taking integer values. Let us briefly 
recall some points of "vocabulary." For a given level a£R,a point t 6 [a, b] 
such that X{t) = a is called "crossing" of the level a. Then Nx(a, [a, b\) 
counts the number of crossings of the level a in the interval [a, b] . Now 
we have to distinguish three different types of crossings (see, e.g., [9]): the 
up-crossings that are points for which X(t) = a and X'(t) > 0, the down- 
crossings that are points for which X(t) = a and X'(t) < and the tangen- 
cies that are points for which X(t) = a and X'(t) = 0. 

Let us also recall that according to Rolle's theorem, whatever the level a 

is, 

Nx(a,[a,b])<N X '(P,[a,b]) + l a.s. 

Note that when there are no tangencies of X' for the level 0, then Nx> (0, [a, b]) 
is the number of local extrema for X, which corresponds to the sum of the 
number of local minima (up zero-crossings of X') and of the number of local 
maxima (down zero-crossings of X'). 
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Dealing with random processes, one may be more interested in the mean 
number of crossings. We will denote by Cx(ot,[a,b]) the mean number of 
crossings of the level a by the process X in [a, b] , 

(3) Cx(a,[a,b])=R(N x (a,[a,b]))=E(#{t€ [a,b] such that X(t) = a}). 

Let us emphasize that this function is no more with integer values and can 
be continuous with respect to a. When, moreover, X is a stationary pro- 
cess, by the additivity of means, we get Cx(a, [a,b]) = (b — a)Cx{a, [0, 1]) 
for a£l. In this case Cx(oc, [0, 1]) corresponds to the mean number of 
crossings of the level a per unit length. Let us also recall that when X 
is a strictly stationary ergodic process, the ergodic theorem states that 
(2T)- 1 iV x (a,[-r,T])^ T ^ +oo Cx(a,[0,l]) a.s. (see [9], e.g.). 

2.1. A Fourier approach for the mean number of crossings function. One 
way to obtain results on crossings for almost every level a is to use the well- 
known co- area formula which is, in fact, valid in the more general framework 
of bounded variations functions (see, e.g., [12]). When X is an almost surely 
continuously differentiable process on [a, b] , for any bounded and continuous 
function /i on 1, we have 

(4) f h(a)N x (a,[a,b])da = [ h(X \t))\X' (t)\ dt a.s. 

JR J a 

In particular, when h = l, this shows that a i— > Nx(a, [a,b]) is integrable on 
K and JmNx(ct, [a,b]) da = J a \X'(t)\ dt is the total variation of X on [a, b\. 
Moreover, taking the expected values we get by Fubini's theorem that 

C x (a,[a,b])da= [ E(\X'(t)\)dt. 



Therefore, when the total variation of X on [a, b] has finite expectation, 
the function a \— > Cx(ot, [a, b]) is integrable on R. This is the case when 
X is also mean square continuously differentiable since then the function 
1 1— > M(\X'(t)\) is continuous on [a, b\. Let us emphasize that this implies, 
in particular, that Cx{ol, [a,b]) < +oo for almost every level a G M but one 
cannot conclude for a fixed given level. However, it allows us to use Fubini's 
theorem such that, taking expectation in (4), for any bounded continuous 
function h, 

(5) / h(a)C x (a,[a,b])da= [ E(h(X(t))\X' (t)\) dt. 



In the following theorem we obtain a closed formula for the Fourier transform 
of the mean number of crossings function, which only involves characteristic 
functions of the process. This can be helpful when considering shot noise 
processes whose characteristic functions are well known. 
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Theorem 1. Let a,b 6l with a <b. Let X be an almost surely and 
mean square continuously differentiable process on [a,b]. Then a h- > Cx (a, 
[a, b]) € L 1 (M) and its Fourier transform Cx(u, [a,b]) is given by 

rb 

(6) C x {u,[a,b})= / E(e iltXW |X'(i)|)dt. 



Moreover, if ipt denotes the joint characteristic function of (X(t), X'(t)), 
then Cx (u, [a, b]) can be computed by 



7T 



1 r b f+°° i /dipt dipt 

C x {u,[a,b}) = / / - — — (u,v) - — — (u,-v) ) dvdt 



1 f^~°° \ 



dv dv 
(ipt(u,v) +ip t (u,-v) -2ip t (u,0))dvdt. 

^ J a JO V 2 

Proof. Choosing in equation (5) h of the form h{x) = exp(iux) for any 
u real, shows that C^(u, [a,b]) = f*E(e iuX W\X'(t)\) dt. Let us now identify 
the right-hand term. Let nt(dx,dy) denote the law of (X (t) , X' (t)) . Then 
the joint characteristic function ip t (u,v) of (X(t),X'(t)) is 

ipt(u,v)=E(exp(iuX(t)+ivX'(t)))= [ e iux+ivy m(dx,dy). 

Since the random vector (X(t), X'(t)) has moments of order two, then tpt is 
twice continuously differentiable on ]R 2 . Now, let us consider the integral 



[ A i(dip t( , dipt, A, 



A r ^y e iux+ivy _ ^y^iux-ivy 

/ nt(dx,dy)dv 

v=0 Jx,y£R 2 v 

-2 f f ye-^l,t(dx,dy)dv 
Jv=0 Jm. 2 v 

ye iux r AVs _^i dv ^ dXjdy y 

Jv=0 V 

The order of integration has been reversed thanks to Fubini's theorem 



x 



\\ye lux sin{ ^ v > | < y 2 which is integrable on [0,A] x R 2 with respect to dv 
Ht(dx,dy), since X'(t) is a second order random variable]. As A goes to 
+oo, then f^ Q sm \"> dv goes to | sign(y), and moreover, for all A, x and y, 
we have \ye lux f v J! sm ^ dv\ < 3\y\, thus by Lebesgue's dominated conver- 
gence theorem, the limit of — -La exists as A goes to infinity and its value 
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is 



1 f A l(d^ t frlH A , 

A->+oo 7T Jq V \ OV OV J 

\y\e iux ^(dx,dy)=E(e mX ^\X'(t)\). 



The second expression in the proposition is simply obtained by integration 
by parts in the above formula. □ 

The last expression considerably simplifies when X is a stationary Gaus- 
sian process almost surely and mean square continuously differentiable on 
R. By independence of X(t) and X'(t) we get tpt(u,v) = <frx(u)<frx'( v ) where 
4>Xi respectively, (ftx', denotes the characteristic function of X(t), respec- 
tively, X'(t) (independent of t by stationarity). Then, the Fourier transform 
of the mean number of crossings function is given by 

C x (u,[a,b]) = -- — -<t>x{u) [ -^^-(v)dv. 

By the inverse Fourier transform we get a weak Rice's formula 

(7) C x (a, [a, &]) = — ( — ) ^ e ^~nx(o))? /2m for ^ Q £ R 
7T \m J 

where m,Q = Var(X(t)) and mi = Var(X'(t)). Let us quote that in fact Rice's 
formula holds for all level a£l and as soon as X is a.s. continuous (see [2], 
Exercise 3.2) in the sense that Cx(&, [o> b]) = +oo if ni2 = +oo. 

However, in general, the knowledge of Cx{u, [«,^]) only allows us to get 
almost everywhere results on Cx(ce, [a, b]) itself, which can still be used in 
practice as explained in [25]. 

2.2. Mean number of crossings for a given level. One way to derive re- 
sults on Cx(ce,[a,b]) for a given level a is to use Kac's counting formula 
(see [2], Lemma 3.1), which we recall now. When X is almost surely contin- 
uously differentiable on [a, b] such that for aGl 

F(3t G [a, b) s.t. X(t) = a and X'(t) = 0) = and 

(8) 

F(X(a) = a)= F(X(b) = a) = 0, 

then, 

1 f b 

(9) N x (a,[a,b]) = \im — J t ]x{t y a]<s \X'(t)\dt a.s. 

The first part of assumption (8) means that the number of tangencies for the 
level q is almost surely. The following proposition gives a simple criterion 
to check this. 
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Proposition 1. Let a, b G R with a<b. Let X be a real valued random 
process almost surely C 2 on [a, b] . Let us assume that there exists cj) £ L 1 (R) 
and c > such that 

Vt€[a,6] \E(e iuX{t) )\ <af>(u). 

Then, 

VaGR P(3te [a,b],X(t) = a and =0) = 0. 

Proof. Let M > and let denote Am the event corresponding to 

max \X"(t)\ < 2M 
te[o,6] 

such that P(3t G [a, b], X(t) = a, X'(t) = 0) = lim M ->+oo G [a,6],X(t) = 
a, = 0, Am)- Let us assume that there exists t G [a, b] such that X(t) = 
a and X'(t) = 0. Then for any n G N there exists £;„ G [2 n a, 2 n b] n Z such 
that |t — 2 _n A; n | < 2 _n and, by the first order Taylor formula, 

(10) \X{2~ n k n ) -a\< 2~ 2n M. 

Therefore, let us denote 

B n = |J {\X(2~ n k n )-a\<2~ 2n M}. 

k n £[2 n a,2 n b]nZ 

Since (B n n j4M)neN is a decreasing sequence we get 

P(3£G [a,b];X(t) = a,X'(t) = 0,A M ) < lim P(B n nA M ). 

n— >+oo 

But, according to assumption, for any n G N the random variable X(2~ n k n ) 
admits a uniformly bounded density function. Therefore, there exists c' > 
such that 

P(|X(2- n A: n ) -a\< 2~ 2n M) < c'2~ 2n M. 
Hence, F(B n n A A/ ) < (6 - a + l)c'2~ n M, which yields the result. □ 

Now taking expectation in (9) gives an upper bound on Cx(cn,[a,b]), 
according to Fatou's lemma, 

C x (a,[a,b])<hminfiy E(l m) _ a \ <s \X' (t)\) dt. 

This upper bound is not very tractable without assumptions on the existence 
of a bounded joint density for the law of (X (t) , X' (t)) . As far as shot noise 
processes are concerned, one can exploit the infinite divisibility property by 
considering the mean number of crossings function of the sum of indepen- 
dent processes. The next proposition gives an upper bound in this setting. 
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Another application of this proposition will be seen in Section 5 where we 
will decompose a shot noise process into the sum of two independent pro- 
cesses (for which crossings are easy to compute) by partitioning the set of 
points of the Poisson process. 

Proposition 2 (Crossings of a sum of independent processes). Leta,b£ 
R with a <b. Let n>2 and Xj be independent real-valued processes almost 
surely and mean square two times continuously differentiable on [a, b] for 
1 < j <n. Assume that there exist constants Cj and probability measures 
d/ij on R such that if dPx j (t) denotes the law of Xj (t) , then 

W G [a, b] dP Xj (t) < c j d Vj f or 1 < j < n - 

Let X be the process obtained by X = Y^=i Xj and assume that X satisfies 
(8) for a G R. Then 

(11) C x (a, [a,b]) < E(]I c M) + 

Moreover, in the case where all the Xj are stationary on R, 

n 

C x (a,[a,b])<^2C x ,(0, [a,b])- 
i=i 

Proof. We first need an elementary result. Let / be a C 1 function on 
[a, b], then for all 5 > 0, and for all x G R, we have 

(12) Y5 f 1 l/(*)-l<5l/'(*)l dt ^ N f'(°> [«' 6 D + L 

This result (that can be found as an exercise at the end of Chapter 3 of [2]) 
can be proved this way: let a\ < ■ ■ ■ < a n denote the points at which f'(t) = 
in [a,b]. On each interval [a,ai], [01,02], . . . , [a n ,6], / is monotonic and 
thus f^ l+1 ^-\f(t)-x\<s\f'{t)\ < 25. Summing up these integrals, we have the 
announced result. 

For the process X , since it satisfies the conditions of Kac's formula (8), 
by (9) and Fatou's lemma, 

Cx(a,[a,6])<hminfiy E(t\X(t)-a\< S \X'(t)\) dt. 
Now, for each 5 > 0, we have 

n 

\X'(t)\) < ^E(l |Xl(4)+ ... +Xn(t) _ a| <,|^(i)|). 
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Then, thanks to the independence of X\ , . . . , X n and to the bound on the 
laws of Xj(t), we get 

6 

^'( 1 \X 1 (t)+-+X n {t)-a\<8\X' 1 {t)\)dt 

/ E ( 1 \X 1 (t)+x 2 +-+x n -a\<s\X' 1 (t)\\ 

X 2 (t)=x 2 ,...,X n (t) = x n ) dP X2{t) (x 2 ), dP Xn (t){x n ) dt 

- J R -if ^hxi(t)+x2+---+x n - a \<s\x[{t)\^ 

Now, (12) holds almost surely for X±, taking expectation we get 
26 

Using the fact the d/ij are probability measures we get 

rb 



1 J m\X 1 {t)+x 2 +- + x n -a\<s\Am dt < C x , (0, [a, b}) + 1. 



— / E(t\Xx(t)+-+Xn(t)-a\<S 

\X[(t)\)dt<(j[c^(C x[ (0,[a,b]) + l). 

3— 2 

We obtain similar bounds for the other terms. Since this holds for all 5 > 0, 
we have the bound (11) on the expected number of crossings of the level a 
by the process X. 

When the Xj are stationary, things become simpler; we can take Cj = 1 
for all 1 < j < n, and also by stationarity we have that for all p > 1 integer 
Cxipt, [a, b + p{b — a)]) = (p + l)Cx(ct, [a, b]). Now, using (11) for all p, then 
dividing by (p + 1), we have that for allp, Cx(ot, [a, b]) < YH=i Cx'. (0) [°> &]) + 
^j-. Finally, letting p go to infinity, we have the result. □ 

As previously seen, taking the expectation in Kac's formula only allows 
us to get an upper bound for C x . However, under stronger assumptions 
(see [18], Theorem 2), one can justify the interversion of the limit and the 
expectation. In particular, one has to assume that (X(t),X'(t)) admits a 
density pt continuous in a neighborhood of {a} x R. Rice's formula states 
that 

(13) C x (a,[a,b])= / / \z\pt(a, z) dz dt < +oo, 

Ja JR 

such that, under appropriate assumptions, one can prove that the mean 
number of crossings function a i— >■ Cx{ol, [a, b]) is continuous on R. 
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3. Crossings of smooth shot noise processes. Prom now on, we focus on a 
shot noise process X given by the formal sum (1), which can also be written 
as the stochastic integral 

(14) X(t)= [ zg{t- s)N{ds,dz), 

JRxR 

where N is a Poisson random measure of intensity \v{ds)F(dz), with F 
the law of the impulse f3 (see [17], Chapter 10, e.g.). We focus in this pa- 
per on stationary shot noise processes for which v(ds) = ds is the Lebesgue 
measure. Such processes are obtained as the almost sure limit of truncated 
shot noise processes denned for UT(ds) = l[_ T ^(s) ds, as T tends to in- 
finity. Therefore, from now on and in all the paper, the measure v{ds) is 
the Lebesgue measure ds or the measure it(gLs). Then, assuming that the 
random impulse /3 is an integrable random variable of L 1 ^) and that the 
kernel function g is an integrable function of L 1 (M), it is enough to ensure 
the almost sure convergence of the infinite sum (see also Campbell's theorem 
and [15]). When, moreover, (3 G L 2 (f2) and g G L 2 (R), the process X defines 
a second order process. 

3.1. Regularity and Fourier transform of the mean number of crossings 
function. Under further regularity assumptions on the kernel function we 
obtain the following sample paths regularity for the shot noise process itself. 

Proposition 3. Let ft G L 2 (0). Letg GC 2 (R) such thatg,g',g" G L 1 (R). 
Then X is almost surely and mean square continuously differentiable on R 
with 

x'{t) = Y J Pi9'{t-T i ) yt G R. 

i 

Proof. Let A > and remark that for any s G R and \t\ < A, since 
geC^R), 

\9(t-s)\= / g'(u- s)du + g(-s) < \g'(u - s)\ du + \g(-s)\, 
Jo J -A 

such that by Fubini's theorem, since g,g' G L X (R), 

/ sup \g(t - s)\ ds < 2A \g'(s)\ds+ / \g(s)\ ds < +oo. 
JR te[—A,A] JR JR 

Therefore, since (3 G L 1 (0), the series A su Pte[-A,A] ~ r «)l converges 
almost surely which means that Y^i ~ T i) converges uniformly on [—A, A] 
almost surely. This implies that the sample paths of X are almost surely con- 
tinuous on R. Similarly, since g' G C X (R) and g',g" G L X (R), almost surely 
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the series X^A</(- — Tj) converges uniformly on [—A, A] and therefore X is 
continuously differentiable on [—.A, A] with X'(t) = Yli PiQ* ~~ r «) f° r au 
i € [—A, A]. Note that the same holds true on [—A + n, A + n] for any 
n € Z, which concludes for the almost sure continuous differentiability on 

K = U«ezM + n.^ + n]- 

Now, let us be concerned with the mean square continuous differentiabil- 
ity. First, g,g' £ L 1 (R) implies that g £ L°°(R) n L 1 (R) C L 2 (R) such that 
X is a second order process since /3 £ L 2 (Q). Its covariance function is given 
by S{t,t') = Cav(X(t),X(t')) = XE{p 2 ) f R g(t - s)g{t' - s)u(ds). Similarly, 
we also have that g' £ L 2 (R) n L°°(R) and X' is a second order process. 
According to [1], Theorem 2.2.2, it is sufficient to remark that assump- 
tions on g ensure that exists and is finite at any point (t, t) £ R 2 with 
$gp(t,t) = AE(/3 2 ) J R g'{t - s)g'{t - s)is(ds). Therefore, for all t G R, the 
limit limfe^o x ( t + h y~ x ( t ) ex ists in L 2 {Q) and is equal to X'(t) by unicity. 
Moreover, the covariance function of X' is given by (t,t') h-> AE(/3 2 ) J R g'(t — 
s)g'{t' - s)is(ds). □ 

Iterating this result one can obtain higher order smoothness properties. In 
particular, it is straightforward to obtain the following result for Gaussian 
kernels. 

Example (Gaussian kernel). Let j3 £ L 2 (J7), g(t) = gi(t) = -h= exp(-i 2 /2) 



and X given by (1). Then, the process X is almost surely and mean square 
smooth on R. Moreover, for any n £ N, 

vt e R = J] /3^ n) (t - r f ) = £ Pi(-i) n H n (t - n) gi (t - Ti ), 

i i 

where H n is the Hermite polynomial of order n. 

From now on, in order to work with almost sure and mean square contin- 
uously differentiable process, we make the following assumption: 

(A) g£C 2 (R) with g,g',g" £L 1 (R). 

Therefore, choosing /3 £ L 2 (Q), the shot noise process X satisfies the as- 
sumptions of Theorem 1 such that the Fourier transform of its mean number 
of crossings function can be written with respect to ipti the joint character- 
istic function of (X(t), X'(t)), given by (see [17], Lemma 10.2, e.g.) 

Vu,v £ R Mu,v) =E(e iuXit)+vX ' (t) ) 

(15) / , \ 

= exp / [ e iz{ug(t-s)+vg'(t-s)) _ l]\ v ^s)F{dz) 
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In order to get stronger results on the mean number of crossings function 
we first have to investigate the existence of a density when considering a 
shot noise process X, or more precisely, a shot noise vector- valued process 
(X, X'). Then we consider an Revalued shot noise process given on M by 

(16) Y(t)=J2PMt-Ti), 

i 

where h : M. i— > M. d is a given (deterministic) measurable vectorial function in 
L 1 (1R). In this setting one can recover X given by (1) with d = 1 and h = g, or 
recover (X, X') (if it exists) with d = 2 and h = (g, g'). It will be particularly 
helpful to see Y as the almost sure limit of a truncated shot noise process 
Yt defined for i>T(ds) = 1[_t,t]( s ) ds, as T > tends to infinity. Therefore, 
from now on and in all the paper, we use the following notation. 

Notation. For any T > 0, we denote by Yt, respectively, Xt, when 
d = l, the shot noise process given by (16), respectively, (1), obtained for 
VT(ds) = l[_ T T ](s) ds. We simply denote by Y, respectively, X, when d= 1, 
the shot noise process obtained for v the Lebesgue measure. 

3.2. Existence of a density and continuity of the mean number of crossings 
function. Let us remark that for d > 1 and T > 0, the shot noise process 
Yt satisfies 

7T 

a?) y t (-)= y, w-^^Ew-^)' 

\n\<T i=i 

where 

(18) lT = #{i;T i £[-T,T}} 

is a Poisson random variable of parameter Ai^t(R) — 2 AT and {U T 1 ^} are 
i.i.d. with uniform law on [—T,T] independent from 77 1 and Here and 

in the sequel the convention is that Y2i=i = an d, as usual, f '= d ' stands for 
the equality in finite dimensional distributions. 

Moreover, for any M > T, one can write Ym as the sum of two indepen- 
dent processes Yt and Ym — Yt such that the existence of a density for the 
random vector Yr(t) implies the existence of a density for the random vector 
Ym($) and therefore for Y(t). Note also that by stationarity Y(s) will also 
admit a density for any sSM. Such a remark can be used, for instance, to es- 
tablish an integral equation to compute or approximate the density in some 
examples [14, 20, 21]. However, the shot noise process may not have a den- 
sity. For example, when h has compact support, there exists A > such that 
h(s) = for \s\ > A. Then, for any T > A, we get X T (0) = X A (0) = X(0) 
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such that P(X T (0) = 0) = P(X(0) = 0) > P(7a = 0) > 0, which proves that 
Xt(0) and X(0) don't have a density. Such a behavior is extremely linked 
to the number of points of the Poisson process {t{\ that are thrown in the 
interval of study. Therefore, by conditioning we obtain the following crite- 
rion. 



Proposition 4. If there exists m > 1 such that for all T > large 
enough, conditionally on {^t — m }> the random variable Yt(0) admits a 
density, then, conditionally on {-jt > m}, the random variable Yr(0) admits 
a density. Moreover, Y(0) admits a density. 

Proof. Let T > be large enough. First, let us remark that condition- 
ally on {jt = Tn}, ^t(O) = Y^iLiPihiU^). Next, notice that if a random 
vector V in M. d admits a density fy then, for Ut with uniform law on [-T, T] 
and j3 with law F, independent of V, the random vector W = V + fih{UT) 
admits ibGR^^ IrI-t fv( w ~ z h(t)) dtF(dz) for density. Therefore, by 

induction, the assumption implies that Y17=l fiih{U^) has a density, for any 
n>m. This proves that, conditionally on {7^ > m}, the random variable 
Yt(0) admits a density. 

To prove that Y(0) admits a density, we follow the same lines as in [3], 
proof of Proposition A. 2. Let A C M. d be a Borel set with Lebesgue measure 
0, since Yt(0) and Y(0) — Yt(0) are independent 

P(F(0) e A) = P(y T (0) + (y(0) - y T (0)) eA)=[ P(y T (0) e A - 

with /z T the law of y(0) - Vt(0). But for any y G K d , 

/ 7t \ 

P(y T (0) 6^-j)=P J /3^(4 4) ) G ^4 — y J 



\i=l 
-00 / 7t 



P ( E W4° )eA- nr = n P( 7T = n) 

n=0 \i=l / 
m—1 / n \ 



n=0 \i=l / 

since A — y has Lebesgue measure and S^Li A(M4 )) nas a density for 
any n > m. Hence, for any T > large enough, 

P(y(0) EA)< P( 7r < m - 1). 

Letting T ->■ +00 we conclude that P(y(0) £ vl) = such that y(0) admits 
a density. □ 
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Let us emphasize that Yt(0) does not admit a density since F(Yt(0) = 
0) > P(7r = 0) > 0. Let us also mention that Breton [6] gives a similar as- 
sumption for real-valued shot noise series in his Proposition 2.1. In particu- 
lar, his Corollary 2.1 can be adapted in our vector- valued setting. 

Corollary 1. Let /i:Rh >M. d be an integrable function and (3 = 1 a.s. 

Let us define hd-^ d >-> M rf by hd(x) = h{x{)-\ \-h(xd), for x = (x±, . . . , xj) S 

M . If the hd image measure of the d-dimensional Lebesgue measure is ab- 
solutely continuous with respect to the d-dimensional Lebesgue measure then 
the random vector Y(0), given by (16), admits a density. 

Proof. Let A C W 1 a Borel set with Lebesgue measure then the as- 
sumptions ensure that j Rd th d (x)eA dx = 0. Therefore, for any T > 0, using 
the notation of Proposition 4, 

P h{uP) E A) = ^ 1^ t hd(x)eA dx = 0. 

Hence, Yli=i h(ufp) admits a density and Proposition 4 gives the conclusion. 

□ 



Example (Gaussian kernel). Let g(t) = -j= exp(— 1 2 /2), (3 = 1 a.s. and 

X given by (1). Let us consider h = (g,g') and h,2 ■ (x±,X2) 6t 2 4 h(xi) + 
h{x2)- The Jacobian of hi is 

J(h 2 )(x 1 ,x 2 ) = — (1 +xix 2 )(xi - x 2 )exp(-(xl + x|)/2). 

Hence, the h 2 image measure of the 2-dimensional Lebesgue measure is 
absolutely continuous with respect to the 2-dimensional Lebesgue measure. 
Then, for any t E M, the law of the random vector (X(t),X'(t)) is absolutely 
continuous with respect to the Lebesgue measure. Note that, in particular, 
this implies the existence of a density for X(t). However, this density is not 
bounded (and therefore not continuous) in a neighborhood of as proved in 
the following proposition. 

Proposition 5. Let us assume for sake of simplicity that (3 = 1 a.s. 
and let g denote the kernel function of the shot noise process. Then: 

1. If g is such that there exist a > 1 and A > such that\/\s\ > A, \g(s)\ < 
e~' s ' a , then 3eo > such that V0 < e < e$, 

P(|X(£)| < e) > \e~ 2XT " where T £ is defined by T £ = (-loge) 1/Q . 
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2. If g is such that there exists A > such that V|s| > A, \g(s)\ < e - '"' 
and if A < l/4 ; t/ien 3eo > such that VO < e < eq, 

* » — 2AT e 



P(|X(i)| < e) > 1-7 tttt e^ AJe w/iere T e is defined by T e = -loge. 

27ms implies in both cases that P(|X(t)| < e)/e goes to +oo as e goes to ; 
and thus the density of X(t) (if it exists) is not bounded in a neighborhood 
0/0. 

Proof. We start with the first case. Let e > and let T £ = (- loge) 1 /". 
Assume that e is small enough to have T e > A. We have by definition 

X(t)^X(0) = Ei9(ri)- If we denote X Te (0)=E| Tl |<T £ 5(r 4 ) and 7^(0) = 
S|tj|>t. 5'( r i)i then Xt s (0) and -Rt £ (0) are independent and X(0) = Xt £ (0) + 
R T J(0). e We also have P(|X(0)| < e) > P(|X Te (0)| = and |i2 Te (0)| < e) = 
P(|X Te (0)| =0)xP(|i2r E (0)| <e).Now, on the one hand, we have P(|X Te (0)| = 
0) > P (there are no n in [-T £ ,T £ ]) = e _2AT£ . On the other hand, the first 
moments of the random variable Rt e (0) are given by E(i?r e (0)) = 
^ f\s\>T 9( s )^ s an< ^ Var(i?T e (0)) = Xf,,™ T g 2 (s)ds. Now, we use the fol- 
lowing inequality on the tail of J e _s " : 



+oo r+oo 

T 



VT>0 e~ Ta = as a - l e~ sa ds>aT a - 1 I e~ s ~ ds. 



Thus, we obtain bounds for the tail of J g and of J g 2 



+oo -T a r+oo P -2T° 

e~ sa ds<^—- and / (e~ sa ) 2 ds< 



aT"" 1 Jt 2aT a ~ 1 

Back to the moments of Rt e (0), since T £ = (— loge) 1//Q we have 

|E(i*r e (0))|<-^r and Var(^(0)) < 



We can take e small enough in such a way that we can assume that 
|E(i?r e (0))| < e. Then, using Chebyshev's inequality, we have 

n\RrM\ < £ ) = - HRtM) < RtM ~ E(i2r e (0)) < e - E(i? Ts (0))) 
> 1-P(\RtM-®(RtM)\ >£- |E(i2 T .(0))|) 
>1 Var^O)) ^ A 



(e- |E(fi T£ (0))|) 2 " * aT £ a -\l - 2A/aT f f- 1 ) 2 

which is larger than 1/2 for T £ large enough (i.e., for e small enough). 

For the second case, we can make exactly the same computations by 
setting a = 1 and get P(|i? Te (0)| < e) > 1 - A/(l - 2A) 2 , which is > when 
A<l/4. □ 
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Such a feature is particularly bothersome when considering crossings of 
these processes since most of known results are based on the existence of 
a bounded density for each marginal of the process. However, this is again 
linked to the number of points of the Poisson process {t{} that are thrown 
in the interval of study. By conditioning, the characteristic functions are 
proved to be integrable such that conditional laws have continuous bounded 
densities. The main tool is Proposition 10 (postponed to the Appendix) 
established using the classical stationary phase estimate for oscillatory inte- 
grals (see, e.g., [26]). 

Proposition 6. Let us assume for sake of simplicity that (3 = 1 a.s., 
let T > 0, a < b and assume that g G L 1 (1R) is a function of class C 2 on 
[-T + a, T + b] such that 



m= min W g' (s) 2 + g" (s) 2 > and 

se[-T+a,T+b] 

(19) 

n = #{sG [-T + a,T + b] s.t. g"(s) = 0} < +oo. 

Then, conditionally on {77 1 > ko} with ko > 3, for all t 6 [a, b] and M > T, 
the law of Xm (t) admits a continuous bounded density. Therefore, for any 
tSl, the law of X(t), conditionally on {jt > ko}, admits a continuous 
bounded density. 

Proof. Actually, we will prove that conditionally on {7^ > £;o}, the 
law of the truncated process X^(t) = Yl\ Ti \<T ~~ r «) admits a continuous 
bounded density for t £ [a,b]. The result will follow, using the fact that for 
M > T, X M (t) = X T (t) + (X M (t) - X T (t)), with X M (t) - X T (t) indepen- 
dent of Xxit). So let us denote Vtfco the characteristic function of Xx(t) 
conditionally on > ko}. Then, for all u £ R, we get 

CcH = E E(e^>|7T = k)F( lT = k) 



fc>fc() 



\ V (— [ T r Wf ~ s) dX r -2AT ( 2 ^ T ) fc 



nir > k 0/ , 

Therefore, 

(20) |V4»|<(2T)- fc ° 



T+t 

e iug{s) ds 

-T+t 



ko 



Hence, using Proposition 10 on [— T + t,T + t] C [— T + a,T + b], one can 
find C a positive constant that depends on T, ko, A, m and no such that for 
any \u\ > 1/m 

\tf M (u)\<C\u\- k °' 2 . 
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Then, since ky > 3, iptk ls integrable on K and thanks to Fourier inverse 
theorem it is the characteristic function of a bounded continuous density. □ 

Using similar ideas we obtain the following result concerning the continu- 
ity of the mean number of crossings function. 



Theorem 2. Assume for sake of simplicity that (3 = 1 a.s. and that g is a 



function of class C 



on . 



satisfying (A). Let T > 0, a < b and assume that for 



all s G [—T + a,T + b], the matrice <&(s) 



wise derivative 3>'(s) = (g(a)?L g( 4 )(sj) are ^ nver tible. Then, conditionally on 
{7t > ^o} with ko > 8, for all M > T, the mean number of crossings function 
a *— > K(Nx M (cx, [a,b])\-jT > &o) ^ s continuous on R. Moreover, 



7 (3) 



and its component- 



is) - 



E(iVx M (a,[a,6])| 7T >fco) 
uniformly on aGl. 



Af- 



E(A^(a,[a,6])| 7 T> fco) 



Proof. The result follows from Rice's formula. To establish it we use [18], 
Theorem 2, and thus we have to check assumptions (i) to (iii) related to joint 
densities. Let t G [a, b] and M > T. We write X M (t) = X T (t) + (X M (t) - 
Xr(t)) with Xm — Xt independent of Xt- We adopt the convention that 
Xoo = X. Let us write for M G [T, +oo] and e small enough 



Vt,e,ko 



, T i T,M 



with ipf^ k the characteristic function of (Xjvf(i), (Xm(1 + s) — XM(t))/s), 
conditionally on > fco}. Note that, — -Xt is independent of 77 1 such 

that ip t ' £ is just the characteristic function of (Xm(^) — -Xt(*)> (O^Af ~~ 
-Xr)(i + £) — (Xm — XT)(t))/e). First we prove that there exists C > such 
that, for all < j < 3, for all M > T and e > small enough, 



(21) 



<C(l + \A* 2 + u 



2 \-(Aso-3)/2 



Let us remark that, since 3', 5" € L^M) by (A), one has g,g' G L°°(R). It 
implies, in particular, that G L 1 (M) n L 2 (M) n L 3 (M) such that the above 
partial derivatives exist. Moreover, by Leibniz formula, for < j < 3, one has 



(22) tt-W 



j\ T M 



1=0 



On the one hand, 



/^,A// \ 



< E 



(X M - X T )(t + e) - (X M - X T )(t) 
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with 



(X M - X T )(t + e)- (Am - X T )(t) 



< Yl \aS-n)\, 

T<\ n \<M 

where g £ (s) = \ f £ g'(s + x)dx is such that g £ G L°° (M)nL 1 (R) with ||# e ||oo < 
H^'lloo and ||g £ ||i < lb' ||i- Then, using the moment formula established in [5], 
one can find c > such that for all < j < 3, with (j — 1) + = max(0, j — 1), 



d 3 ,T,M ( \ 



(23) 



<e(Y Yl Ideit-T^Y) 

< cmax(l, H^lloo) '- 1 ^ max(l, Alb'lUP. 



On the other hand, 

P(7T> MC.fcoM 



k>ko 
k>ko 



where 



X Uu,v) = (2T)- 1 



T+t 



T+t 



Jug(s)+ivg e (s) ^ g 



is the characteristic function of (g(t — UT),g e (t — Ut)), with Ut a uniform 
random variable on [—T,T]. It follows that \xt £ {u,v)\ < 1, so that one can 
find c > such that for all < j < 3, 



7^ e ,fco(«.w) 



< C max(l,|| 5 '|| 00 )^ 1 )+max(l,A||< 7 '|| 1 y 

- P(7r - fc °" j) lxL(M)l^. 



P( 7T > fco) 

This, together with (23) and (22), implies that one can find c> such that 
for all < j < 3, 



(24) 



<9 j 



< C max(l,|| 5 / || 00 )(^ 1 )+max(l,A||< 7 / || i y 



P( 7T > fe ) 



Moreover, let Q £ {s 



) and $' e ( S ) = 1^). Then det$ £ ( S ) 



Voo fl?(«) 

converges to det $(s) as e — > 0, uniformly in s 6 [— T — a, T + 6]. The assump 
tion on $ ensures that one can find £q such that for e < £q, the matrix 3> e (s) 
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is invertible for all s G [-T — a,T + b]. The same holds true for $j(s). Denote 
m = Toam-se[-T-a,T+b],e<e \\^s( s ) > where || • || is the matricial norm 
induced by the Euclidean one. According to Proposition 10 with no = 0, 



V(u, v) G R 2 s.t. \Jv? + v 2 > — 



\ X Uu,v)\ = (2T)- 



in 

T+t 



e iug(s)+ivg E (s) ^ g 

T+t 



< 



24^/2 



V my/u 2 + v 2 

Therefore, one can find a constant c ko > such that, for all < j < 3, 

I I is less than 

c fc0 (2r)- fc0+3 max(l,|b / || oo )^- 1 )+max(l,A|| 5 / || 1 )^ 



x 



P(7T>*ti-j) (1 + ,/-^ r (*o-3)/2 



P(7T > ko) 



Letting e tend to we obtain the same bounds as (21) for ip^ ko the charac- 
teristic function of (XM{t), X' M (t)) conditionally on {^t > &o}- Since ko > 8, 
(21) for j = ensures that iff Ejko G L^M 2 ), respectively, i t A f Q G L X (R 2 ), 
such that, conditionally on {7-^ > ko}, (Xm(£), (^m(* + e) — ^M(i))/&)) re- 
spectively, (X M (t), admits p&Jx.z) = ^ / R2 e--™^ fc >, 

v)dudv, respectively, p^ ko = 4^ J R 2 e~ lxu ~ lzv ipfl. o (u,v) dudv, as density. 
Moreover: 

(i) Pt e ko (x, z) is continuous in (t, x) for each z,e, according to Lebesgue's 
dominated convergence theorem using the fact that Xm is almost surely 
continuous on R. 

(ii) Since Xm is almost surely continuously differentiable on R we clearly 
have for any (it, v) G R 2 , Vfe ko^ u i v ) ~~ ^ ^t k i u ^ v ) as e — » 0. Then by Lebes- 
gue's dominated convergence theorem, using (21) for j = we check that 
P^e k ( x i z ) ~~ ^ Ptk z ) as e — )• 0, uniformly in (t, x) for each z G R. 

(iii) For any z 7^ 0, integrating by parts we get 



= ^273 / . ^ XU ~ lZV T^^,k (u,v)dudv, 



X 



such that by (21) for j = 3, we check that p^ £ k (x,z) < Ch(z) for all £,e 
with /i(z) = (1 + |z| 3 ) -1 satisfying J* K |z|/i(z) dz < +00 and C a positive con- 
stant. 

Therefore, [18], Theorem 2, implies that 

E(N XM (a, [a,b])\7T >k ) = [ f \z\p™ ko (a, z) dzdt, 

X 



which concludes the proof, using pf ko (a,z) = A7T l z s | s2 e iau izv -^s 
tp^ ko (u,v) dudv and (21) for j = 3. □ 
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Note that, despite that we have closed forms, these crossings formulas are 
not very tractable for general shot noise processes. However, as the intensity 
A of the shot noise process X tends to infinity, due to its infinitely divisible 
property and since it is of second order, we obtain, after renormalization, a 
Gaussian process at the limit. It is then natural to hope for the same kind 
of asymptotics for the mean number of crossings function. This behavior is 
studied in detail in the next section. 

4. High intensity and Gaussian field. 

4.1. General feature. It is well known that, as the intensity A of the 
Poisson process goes to infinity, the shot noise process converges to a nor- 
mal process. Precise bounds on the distance between the law of X(t) and 
the normal distribution are given by Papoulis [22]. Moreover, Heinrich and 
Schmidt [15] give conditions of normal convergence for a wide class of shot 
noise processes (not restricted to processes defined on R, nor to Poisson 
processes). In this section we obtain a stronger result for smooth station- 
ary shot noise processes by considering convergence in law in the space of 
continuous functions. In all of this section we continue to assume that X is 
a stationary shot noise process obtained for v the Lebesgue measure on R, 
and we will denote X\ the strictly stationary shot noise process given by (1) 
with intensity A > 0. Let us define the normalized shot noise process 

(25) z x (t) = ^=(x x (t)-E(x x (t))), tern. 

Then, we obtain the following result. 

Proposition 7. Let (3 e L 2 (£l) and g satisfying (A). Then, 




where B is a stationary centered Gaussian process almost surely and mean 
square continuously differentiable, with covariance function 

Cov(B(t),B(t'))= [ g{t-s)g{t' -s)ds. 
Jr 

When, moreover, g" 6 L P (R) for p > 1, the convergence holds in distribu- 
tion on the space of continuous functions on compact sets endowed with the 
topology of the uniform convergence. 

Proof. We begin with the proof of the finite dimensional distributions 
convergence. Let k be an integer with k > 1 and let t\, . . . , t k S R and w\ = 
(ui,vi),...,w k = (u k ,v k ) GR 2 . 
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Let us write 

k 



7=1 V \ j \ j / / 



for g(s) = Yl k j=i{ u j9{tj — s) + Vjg'(tj — s)). Therefore 



bgE(e*^=i yA( **^) = A 
Note that as A — > +00 , 



iz(g(s)/VX) 



X ( e iz(g(s)/VX)_ 1 _ iz 9i^ 

V A 



with, for all A > 0, 

Aexpl 



1 



\/a / Va 

By the dominated convergence theorem, since g £ L 2 



-1) 



2~Y \2 



dsF(dz) 



and /3 £ L 2 (£l), we 



get that, as A — > +00, 

k 



E 



exp ij^ 



exp 



3=1 



^E(/3 2 ) / g(sfds). 



Let us identify the limiting process. Let us recall that X\ is a second order 
process with covariance function given by Cov(X\(t), X\(t')) = XE((3 2 )S(t — 
t') with S(t) = f R g(t — s)g(—s) ds. Hence, one can define B to be a stationary 
Gaussian centered process with (t, t') h-» S(t — t') as covariance function. The 
assumptions on g ensure that the function S is twice differentiate. There- 
fore B is mean square differentiate with B' a stationary Gaussian centered 
process with (t,t') h-» —S"(t — t') = f R g'{t — t' — s)g' '(— s) ds as covariance 
function. Moreover, 

E((B'(t) - B\t')) 2 ) = 2(5" (0) - 5"(t - 0) < 2||^||oo||/||i|< - 

such that by [1], Theorem 3.4.1., the process B' is almost surely continuous 
on R. Therefore, as in [11], page 536, one can check that almost surely B{t) = 
B(0) + J B'(s)ds, such that B is almost surely continuously differ entiable. 
We conclude for the f.d.d. convergence by noticing that 

k 



/ g(s) 2 ds = Var ( V UjB(tj) + VjB'itj] 

Jr \^ 



Let us prove the convergence in distribution on the space of continuous 
functions on compact sets endowed with the topology of the uniform con- 
vergence. It is enough to prove the tightness of the sequence (1a) a according 



22 



H. BIERME AND A. DESOLNEUX 



to [17], Lemma 14.2 and Theorem 14.3. Let t, s E R and remark that for any 
q > 1, on the one hand, 

E((Z A (t) - Z x (t')f) = E((3 2 ) [ (g(t -s)- g(t' - s)) 2 ds 

JR 

<E(/3 2 )|| 5 '||J«/|| 1 |i-i'| 2 -V<z. 

On the other hand, 

E{{Z' x {t) - Z' x {t')f) = E{f) [ ( g '(t-s)-g'(t'-s)) 2 ds 

JR 

<E03 2 )||/||,||/|| 1 |t-*'| 2 - 1 /9. 

Note that, assuming that g" G L^R), it allows us to choose q = p > 1 in the 
second upper bound such that 2 — l/q> 1. Moreover, assumption (A) implies 
that g' G L°°(R) n L X (R) C L P (R) such that one can also choose q = p in the 
first upper bound. Then, (1a) A satisfies a Kolmogorov-Chentsov criterion 
which implies its tightness according to [17], Corollary 14.9. □ 

In particular, when a<b, the functional (/, g) \-t h(f(t))\g(t)\ dt is 
clearly continuous and bounded on C([a,6],R) x C([a, 6],R) for any continu- 
ous bounded function h on R. Then, Proposition 7 implies that 

f\(h(Z x (t))\Z' x (t)\)dt — > [ b E(h(B(t))\B'(t)\)dt. 

By the co-area formula (4), this means the weak convergence of the mean 
number of crossings function, that is, 

c z x (-, [a,b]) ^a^+oo C B (; [a,b]). 

This implies also the pointwise convergence of Fourier transforms. Such a 
result can be compared to the classical central limit theorem. Numerous 
improved results can be obtained under stronger assumptions than the clas- 
sical ones. This is the case, for instance, for the rate of convergence derived 
by the Berry-Esseen theorem or the convergence of the densities. We refer 
to [13], Chapters 15 and 16. Adapting the technical proofs allows us to get 
similar results for crossings in the next section. 

4.2. High intensity: rate of convergence for the mean number of crossings 
function. Let us remark that only E(/3 2 ) appears in the limit field. For 
sake of simplicity we may assume that /3 = 1 a.s. Note that, according to 
Rice's formula [9], as recalled in equation (7), since the limit Gaussian field 
is stationary, Cb{ch, [a, b]) = {b — a)Cs(a, [0, 1]) with 

i / \ V 2 

C B (a, [0,1]) = -( — e -« 2 /2m Va G R, 
7r \m J 
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where ijiq = Vax(B(t)) = f R g(s) 2 ds and = Var(_B'(t)) = f R g'(s) 2 ds. More- 



over, its Fourier transform is given by Cb(u, [0, 1]) = y^^-e m ° u ^ • We ob- 
tain the following rate of convergence, for which the proof is postponed to 
the Appendix. 

Proposition 8. Let (3 = 1 a.s. and let g satisfy (A). There exist three 
constants a±, ai and a 3 (depending only on g and its derivative) such that 

VA > 0, Mu G R such that \u\ < ai\^X then 



C? A (u,[0,l])-y^e-^ 2 / 2 
where tuq = f R g(s) 2 ds and 1712 = J R g'(s) 2 ds. 



a 2 + a 3 \u\ 



Let us emphasize that this implies the uniform convergence of the Fourier 
transform of the mean number of crossings functions on any fixed interval. 
Moreover, taking u = 0, the previous upper bound may be a bit refined such 
that the following corollary is in force. 

Corollary 2. Let (3 = 1 a.s. and let g satisfy (A). The mean total 
variation of the process satisfies 



VAX) 



E(\X' x (t)\) 2m 2 



V\ V TT 

where m,2 = J R g'(s) 2 ds and 771,3 = Jjj lff'( s )| 3 ds. 



147723 
37r7772 1 /A ' 



Under additional assumptions we obtain the following uniform conver- 
gence for the mean number of crossings function. The proof is inspired 
by [13], Theorem 2, page 516, concerning the central limit theorem for den- 
sities. 

Theorem 3. Let (3 = 1 a.s. Let us assume, moreover, that g is a func- 
tion of class C on K satisfying (A) such that for all s £ [—1,2], &(s) = 
(g'{s) g"(s) x , $/ , s _ / g "(s) g^(s)\ invertible 

Let 7a = 4f{i]T\^ G [—1, 1]} with {r\ t i}i the points of a Poisson point pro- 
cess with intensity A > 0. 
Then 

C Zx (a, [0, 1]| 7 a > A) — > C B (a, [0, 1]) = - (^) ^ \-°?l 2 ™* 

A-s-+oo TT \777o / 

uniformly in a G M, 

where niQ = J R g(s) 2 ds and 7712 = J R g' (s) 2 ds. 
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Proof. Let A > 8. Then, according to Theorem 1, Cz x (u, [0, 1] > A) 
and Cb(u, [0, 1]) are integrable such that Cz x (a, [0, 1]| 7a > A) and Cb(q, [0, 1]) 
are bounded continuous functions with, for any a£R, 

\C Zx (a,[0,l]\ lx > X) - C B (a,[0,l])\ 

I \Cz~ x (u,[0,l]hx> X) - CB(u,[0,l})\du. 
to Jr 

Let then 

Czl(u, [0,l])-C? A (u, [0,1]| 7a >A) 

E (e^(°)|^(0)|l 7A<A ) - S?^C^(n, [0,1]). 



P( 7a >A) v 1 AWI ^ P( 7 a>A) 

Note that \Cz x (u, [0, 1])| < E(|Z A (0)|), which is bounded according to Corol- 
lary 2, while by the Cauchy-Schwarz inequality, 

|E( e »^°)|Z A (0)|l 7A<A )| < E(Z A (0) 2 ) 1 / 2 P( 7A < A) 1 / 2 , 

with E(Z A (0) 2 ) = Var(Z A (0)) < max(l, ||p'|| 0O )||s'||i. Therefore, one can find 
c\ > such that 

— — Fh x < X) 1/2 

\C Zx (u, [0, 1]) - C Zx (u, [0, 1]| 7 a > A) | < Cl ^ > A) • 

According to Markov's inequality, 

P( 7A < A)=P( e -M2)7X > e -ln(2)A) < E ( e -ln(2)( 7 A-A)j =exp (_(l _ l n (2))A). 

Choosing A large enough such that, in particular, ^p^^j < 77^ according 
to Proposition 8 one can find c 2 such that for all \u\ < A 1 / 8 , 

|C? A (u, [0, 1]| 7 a) - Cb(u, [0, 1])| < c 2 A~ 3 / 8 . 
Thus we may conclude that 

[0, 1]|7A > A) - Cb(u, [0, 1])| du — > 0. 

| u | <A i/8 A^+oo 

Now, let us be concerned with the remaining integral for |u| > A 1 / 8 . Accord- 
ing to Theorem 1, 

Cz A (n, [0, 1]| 7 a > A) = Cx A f [0, 1]| 7 A > A 

with C? A (^, [0, 1]| 7 A > A) = J 1 E(e J W^)^( i )|X A (t)||7 A > X)dt and 
E(e^/ A/X ^W|X A (i)|| 7A >A) 



1 



1 ( dipt,x ( u v \ dip^x ( u v 



+oo 

dv, 



ir J v \ dv Vv 7 ^' VxJ dv Wa' v 7 ^ 
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where ipt,\ is the characteristic function of (X\(t), X' x (t)) conditionally on 
{7a > A}. Integrating by parts we obtain 

- 1 1 ( dvb t \ ( u v 



v \ dv 



dipt,\ ( u v 
dv V^A' VX 



dv 



u v \ 



dv. 



Then, according to (24), one can find a positive constant C3 > such that 

\E{e< u ^ x *®\X' x (t)\\ix> A) I 
v2 P( 7A >A-2) 



P(7A > A) 



Xt 



u v 

IVl'Tx 



A-2 



-^=1 MkDI^o^lwl^l + W l ^\v\>i ) dv 



where Xt{u, v) = | J^^L e m f( s )+ iv s'( s ) <fs is the characteristic function of (g(t - 
U),g'(t — U)), with U a uniform random variable on [—1,1]. Then, 

|C^>,[0,1]| 7 a>A)-Cb(u, [0,l})\ du 



H>avs 

<f |^(u,[0,l]|7A>A)|d«+ / |^(«,[0,l])|d« 

J|«|>AV8 J|u|>AV8 

= /i(A) + J 2 (A). 

Now, for 6 G [0,27r], let us consider the random variable V^g = cos(9)g(t 
U)+sm(6)g'(t-U) such that for any r > 0, Xt(rcos(6),rsm(6)) =E(e iH/ *.») 
(pt,e( r )- By a change of variables in polar coordinates, since A > 1, we get 

A-2 



/i(A)<c 4 (A) 



-00 p2n 



A 1 / 8 JO 



<Pt,e 



r(|ln(r|sin(0)|)| +l)d0dr, 



with c 4 (A) = c 3 A 3 / 2 P {7 ( ^>^ 2) . Since det$(s) ^ for any s G [-1 + i, 1 + t], 
we have the following property (see [13], page 516): there exists 5 > such 
that 

\tf>tfi(r)\ < e -M*)/ 4 ) r2 Vr G (0,5], V0 G [0,2tt] and 
77= sup \ip t) e(r)\<l, 

r>S,e&[0,2n] 

with = minegro^Tr] Var(T4,e) > 0. Note also that according to Proposi- 
tion 10, |(^ ti e(r)| < 24y / ^|r~ 1 / 2 for any r > m with m = min s6 [_ 12 ] ll^ > ( s ) _1 II _1 > 
which may be assumed to be larger than 5. Then, for A large enough such 
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that A 1 / 8 G (e,<h/A), 

/ e- (K(t)/8)Al/4 rln(r)(ir+ / r/ A " 2 r ln(r) 



5 

f +00 



+ ( 24 v^j I^" A " 7r " 3/2l,,(r) * / 

with C5(A) = C4(A)(/ Q 27r (2 + | ln(| sin(0)|)|) d6). This enables us to conclude 
that Ii(A) — ^-j.+oo 0. This concludes the proof since clearly /2(A) — ^,\_s > _|_ 00 
0. □ 

Notice that to obtain the convergence in Theorem 3 without the condi- 
tioning on {7^ > A} (which is an event of probability going to 1 exponentially 
fast as A goes to infinity), one simply needs to have an upper-bound poly- 
nomial in A on the second moment of the number of crossings Nz x (ot, [0, 1]). 

5. The Gaussian kernel. In this section we will be interested in a real 
application of shot noise processes in physics. Indeed, each time a physical 
model is given by sources that produce each a potential in such a way that 
the global potential at a point is the sum of all the individual potentials, 
then this can be modeled as a shot noise process. In particular, we will be 
interested here in the temperature produced by sources of heat. Assuming 
that the sources are randomly placed as a Poisson point process of intensity 
A on the real line R, then the temperature after a time a 2 on the line is 
given by the following shot noise process X\y. 

t G M m. X Ka {t) = -4=e-(^) 2 /^ 2 , 
cry 2tt 



where the {rj} are the points of a Poisson process of intensity A > on R. In 
the following, we will denote by g a the Gaussian kernel of width a defined 
for all t G R by 

9.® = -^—^'^ 



crV2vr 

We will be interested in the crossings of X\^ a because they provide infor- 
mation on the way the temperature is distributed on the line. The number 
of local extrema of X\ a is also interesting for practical applications since it 
measures the way the temperature fluctuates on the line. In a first part, we 
will be interested in the crossings of X\, a when A and a are fixed, and then, 
in a second part, we will study how the number of crossings evolves when 
these two parameters change. From the point of view of applications, this 
amounts to describing the fluctuations of the temperature on the line when 
the time (recall that a 2 represents the time) increases, or when the number 
of sources changes. 
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5.1. Crossings and local extrema of X\ t(T Crossings and local extrema of 
X lambda, sigma. We assume in this subsection that A > and a > are 
fixed. Since the Gaussian kernel g a , and its derivatives are smooth functions 
which belong to all L p spaces, many results of the previous sections about 
crossings can be applied here. In particular, we have: 

• the function a \— > Cx x ^(ct, [a, b]) belongs to L 1 (M) (by Theorem 1); 

• for any T > 0, the function a h-> Cx x a (ot, [a, b\\^T > 8) is continuous (by 
Theorem 2), with 7T = #{n £ [-T,T]}. 

This second point comes from the fact that the Gaussian kernel satisfies 
the hypothesis of Theorem 2. Indeed, the derivatives of g a are given by 

Qa{s) = a ^2^ e ~ s2 ^ 2a2 ' ^~"fe -Hfc(f); where the H^s are the Hermite poly- 
nomials (H\(x) = x ; H2(x) = x 2 — 1; H-^(x) = x 3 — 3x and H^{x) = x 4 — 
6x 2 + 3). Thus, using the notation of Theorem 2, we get det$(s) = ^-(^ + 
l)(^=e~ s2 /2- 2 )2 < o and de t$'( s ) = 4(4 + 3 )( i e - s2 / 2ff2 ) 2 < 0. These 
two matrices are thus invertible for all s 6 1R. 

The first point implies that for almost every a£l, the expected num- 
ber of crossings of the level a by X\ )a is finite. We will now prove in the 
following proposition that in fact, for every a € E, Cx XrT (ot, [a,b]) < +oo, 
by considering the zero-crossings of the derivative X' x and using Rolle's 
theorem. 

In the sequel, we will denote by p(A, cr) the mean number of local extrema 
of X\ <a in the interval [0, 1]. It is the mean number of local extrema per unit 
length. 

Proposition 9. We have 

F(3t £ [0, 1] such that X' X (T (t) = and X" tCT {t) = 0) = 0, 

which implies that the local extrema of X\^ a are exactly the points where the 
derivative vanishes; in other words p(\,cr) =K(Nx' (0, [0, 1])). Moreover, 
we have the following bounds: 

V«GM Cx Xiir (a,[0,i\) <p(X,a) < (3A(2 + 2a) + l)e A . 

Proof. For the first part of the proposition, we use Proposition 10 (in 
the Appendix) with the kernel function h = g' a on the interval [— T + 1, T] for 
T > 0. For this function we can compute h'(s) = ^ 3 1 / ^- (— 1 + ^)e~ s2 ^ 2cr2 and 
h"(s) = ^7^(3f - ^)e~ s2 / 2a2 , and thus * n Q = 3 and m(a,T) = 
min s6 [_< rjr+1 ] \J h'(s) 2 + h"(s) 2 > (we do not need to have an exact value 
for it but notice that it is of the order of e _T2 / 2<j2 when T is large). Finally, 
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as in (20), we get that there is a constant c(T, a) which depends continuously 
on a and T such that 

^3 



|E(e |tt — 3)| ^ — = 

(i + v Mr 



with 7r = #{Tj € [— T, T]}. We can now use Proposition 1 and we get that 
for all T> 1, 

P(3t G [0, 1] such that X' Xa (t) = and X'{ a {t) = Q\ lT > 3) = 0. 

Since the events {75" > 3} are an increasing sequence of events such that 
P(7t > 3) goes to 1 as T goes to infinity, we obtain that ¥(3t G [0, 1] such that 
X^(*) = and = 0)=0. 

For the second part of the proposition, the left-hand inequality is simply 
a consequence of Proposition 2 for the process X' x and n = 1. 

To obtain the right-hand inequality [the bound on p(X,a)], we will apply 
Proposition 2 to the process X' x for the crossings of the level on the 
interval [0, 1]. We already know by the first part of the proposition and by 
Corollary 1 that condition (8) for Kac's formula is satisfied by X( . Then 
we write, for all t G [0, 1], 

xi,„( t ) = £^-n) * E z&pl e -(,-n^ 

+ 1 V" ~( t ~ T ») c -(t-r,;) 2 /2<T 2 

Let [resp., ^(t)] denote the first (resp., second) term. We then have 

Since (i-r^) 2 > a 2 for all t G [0, 1] and all n G R\ [-a,l + a], we get Y£{t) > 
on [0,1] and thus iVy/(0, [0, 1]) =0 a.s. Note that when the event #{rj G 
[-a, 1 + a]} = holds, then = Y" 2 such that N x > (0, [0, 1]) < 1. On 
the other hand, let us work conditionally on #{rj G [— a, 1 + <r]} > 1. The 
probability of this event is 1 — e _A ( 1+2<T ). To study the zero-crossings of Y{, 
we first need an elementary lemma. 

Lemma 1. Let n> 1 be an integer. Let P\,...,P n be n real nonzero 
polynomials and let a%, . . . ,a n be n real numbers, then 

{n *\ n 

t£m. such that ^Pi(t)e a ^ = \ < ^deg(Pj) + n - 1. 
i=l J i=l 
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This elementary result can be proved by induction on n. For n = 1, it 
is obviously true. Assume the result holds for n > 1, then we prove it 
for n + 1 in the following way. For t G R, Ya=i Pi(t)e ait = <=^ /(t) := 
P n +i{t) + TJi=i Pi(t)e( ai ~ an + 1 ^ t = 0. Let k denote the degree of P n+1 . Thanks 
to Rolle's theorem, we have that N f (0,R) < N r (0,R) + 1 < iV///(0,R) + 
2 < ■ ■ < iV /(fc+ i)(0,R) + fc + 1. But can be written as = 

X^=l Qi{t)z < " ai ~ an+1 ^ t i where the Qi are polynomials of degree deg(Qj) < 
deg(Pj). Thus by induction Nj-( k +i)(0,R) < ^r=ideg(-Pi) + n — L and then 

Nf(0, R) < ELi deg(Pi) + n - 1 + k + 1 < ^=1 deg(Pi) + n. This proves 
the result for n + 1. 

Thanks to this lemma, we get that iVy/ (0, [0, 1]) < 3#{n G [-a, 1 + cr]} - 1 
such that 

E(iV y/ (0, [0, l])|#{ri G [-a, 1 + a}} > 1) < 3A(1 + 2<x)/(l - e " A ( 1+2<T )) - 1. 

To use Proposition 2, we need to obtain uniform bounds on the laws of 
Y\{t) and of ^(i) when t G [0, 1]. As in the notation of the proposition, we 
will denote these constants by c\ and C2- Let us start with Y\. Let U be a 
random variable following the uniform distribution on [— 1 — a, 1 + a] . For 
t G [0, 1], we can write U as U = i] t Ut + (1 — r] t )Vt, where Ut is uniform on 
[—l — a + t,a + t], Vt is uniform on [— 1 — a, — 1 — a + t] U [cr + t, a + 1] and % 
is an independent Bernoulli random variable with parameter 2+20 m ^ e then 
have g' a (U) = rj t g' a (Ut) + (1 — r)t)g' a (Vt). Thus the law of g' a (U) is the mixture 
of the law of g* a (Ut) and of the one of g' a (Vt), with respective weights 5+2? 
and 1 — • Consequently 

ViG [0,1], VxGM ^(C/ t )N <Y^dP g , AU) {x). 
The law of Yi(t) conditionally on #{Tj G [— a, 1 + cr]} > 1 can be written as 

dP Yl (t)(x) 

Thus, if we write /o = dPg^^, we get 

cWYi(t)0c) 



- 1 _ e -A(i+ar) fc j ll + 2a/ 

fc=i x 



(/o* ••■ */o)(a:) 



1 _ e -A(2+2 ( r) _ 
6 l_ e -A(l + 2.) /o( X )' 
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where fo(x)dx is a probability measure on R. This shows that we can take 

_ A 1 _ e -A(2+2<x) 
C l — e 1 _ e -A(l+2 c r) • 

For Y2(t), we first notice that Y2(t) can be decomposed as the sum of two 
independent random variables in the following way: 

Y 2 (t)= Yl aUt-n) 

Ti6(— oo,-l — o-+t]U[l+<7+t,+oo) 

+ g'At-n). 

ne(-<r-i+t,-tr)u(i+(T,i+(T+t) 

The first random variable in the sum above has a law that does not depend 
on t. For the second random variable, using the same trick as above [i.e., 
decompose here a uniform random variable on the interval ( — 1 — a,— a) U 
(ct,1 + o~) as a mixture with weights 1/2 and 1/2 of two uniform random 
variables: one on (— 1 — a, — 1 — a + t) U (t + a, 1 + a), and the other one on 
the rest] , we obtain that C2 = e A . 

And finally the bound on the expectation of the number of local extrema 

is 



p(A,cr) < ci 



3A(1 + 2(J ) + C9 ^ fl _ e -MWh + e -A(l+2a) 
1 _ e -A(l+2a) ^ C2 J ^ e >^ e 



< e A ^^(3A(l + 2a)) + e x = (3A(2 + 2a) + l)e A . Q 

5.2. Scaling properties. An interesting property of the shot noise process 
with Gaussian kernel is that we have two scale parameters: the intensity A 
of the Poisson point process and the width a of the Gaussian kernel. These 
two parameters are linked in the sense that changing one of them amounts 
to change the other one in an appropriate way. These scaling properties are 
described more precisely in the following lemma. 

Lemma 2. We have the following scaling properties for the process X\ a : 

1. Changing a and X in a proportional way: for all c > 0, 

{X A/C)CCT (i); t G M} f | \ X A , CT ; t G R J. 

2. Increasing the width of the Gaussian kernel: for all o~\ and a2, we have 

{A aV ^ G R} a = S - {(X x>ai * g a2 )(t);t G R}. 

3. Increasing the intensity of the Poisson process: for all c > 0, we have 
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4. The mean number p(X,o~) of local extrema of X\ >a per unit length sat- 
isfies 

Vc> cp(A, ccj) = p(c\, a). 

Proof. For the first property, let {t^} be a Poisson point process of 
intensity A/c on the line. Then 

*W<> = £ -W- (, -" ,V2cV = ; E * f ; - - 



cay 2vr 

Since the points {t«/c} are now the points of a Poisson process on intensity 
A on the line, we obtain the first scaling property. The second property 
comes simply from the fact that if g ai and g a2 are two Gaussian kernels of 
respective width a\ and o"2, then their convolution is the Gaussian kernel of 
width yj o~\ + o\ . The third property is just a consequence of combining the 
first and second properties. 

For the fourth property, we first compute 

1 y^ -(*-T j) 



^\,co\ L ) - f3. n 2 e 



1 -(t/c-Tj/c) c _ (t/c _ n/c) 2 /2a 2 



c 2 o-\/2k o~ 

where the {r,} are the points of a Poisson point process of intensity A on R. 
Then, since the {Tj/c} are now the points of a Poisson point process of 
intensity cA on R, we have that the expected number of points t £ [0, c] 
such that X' Xca (t) = [which, by definition, equals cp(X,ca)], also equals 

the expected number of points t G [0,1] such that X' cXa (t) = [which is 
p(cX,a)]. □ 

To study how p(X, a) varies when A and a vary, we first can use the result 
on high intensity and convergence to the crossings of a Gaussian process 
obtained in Theorem 3. Indeed, if the second moment of A^' (0) is bounded 
by a polynomial in A, then we will get 

And thanks to the scaling properties, this also will imply that p(A, a) is 
equivalent to ^y^f as a goes to +oo. These two facts have been empirically 
checked and are illustrated on Figure 1. Now, notice that we can also observe 
on the left-hand figure another regime when A is small. Indeed, p(X, a) seems 
to be almost linear for small values of A. Notice also on the right-hand figure 
that p(X,o~) seems to be a decreasing function of a (this then indicates that, 
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0.5 1 
Intensity lambda of the Poisson Point Process 



1 1.5 2 2.5 3 3.5 
Width a of the Gaussian kernel 



Fig. 1. On the left: empirical mean number of local extrema of X\^ a per unit length as 
a function of A (here a = 1 and we have taken the mean value from 50 samples on the 

interval [—100,100]). The horizontal dashed line is the constant ^\J\ and the dotted line 
is the map A i— ► 2A. On the right: empirical mean number of local extrema of X\ ttJ per unit 
length as a function of a (here A = 1 and we have taken the mean value from 10 samples 
on the interval [—100,100]). 



as time goes by, the temperature on the line fluctuates less and less). The 
study of these two facts is the aim of the next section. 

5.3. Heat equation and local extrema. In this subsection we assume first 
that A > is fixed. As we already mentioned it in the introduction of Sec- 
tion 5, one of the main features of the shot noise process X\ a is that it can 
be seen in a dynamic way, which means that we can study how it evolves 
as the width a of the Gaussian kernel changes and consider it as a random 
field indexed by the variable (cr,t). Then, the main tool is the heat equation 
which is satisfied by the Gaussian kernel 



(26) 



Vcr > 0, Vt G 



~^~(^) = u 9a(t) an d a ^ so consequently = a 9^\^)- 



Since the Gaussian kernel g a is a very smooth function, both in a > and 
t G R, by the same type of proof as the ones in Proposition 3, we have that 
(o~,t) !->■ Xx >( r(t) is almost surely and mean square smooth on (0,+oo) x R 
with 



Xx -^-(t) = ^ -j E -( t - T i)= Q-X'{ a {t) and also 



da 



(27) 



da 



X ^t)=aX^(t). 



We will see in the following that this equation will be of great interest to 
study the crossings of X\ a . 
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The convolution of a real function denned on R with a Gaussian ker- 
nel of increasing width a (which amounts to apply the heat equation) is 
a very common smoothing technique in signal processing. One of its main 
properties is generally formulated by the wide-spread idea that "Gaussian 
convolution on R cannot create new extrema" (and it is in some sense the 
only kernel that has this property; see [27]). This has been studied (together 
with its extension in higher dimension) for applications in image process- 
ing by Lindeberg [19], and also by other authors (e.g., to study mixtures of 
Gaussian distributions as in [8] and [7]). However, in most cases, the cor- 
rect mathematical framework for the validity of this property is not exactly 
stated. Thus we start here with a lemma giving the conditions under which 
one can obtain properties for the zero-crossings of a function solution of 
the heat equation. The result, which proof is postponed to the Appendix, is 
stated under a general form for a function h in the two variables a and t. But 
we have to keep in mind that we will want to apply this to h(a, t) = X' x (t) 
to follow the local extrema of the shot noise process when a evolves. 

Lemma 3. Let ao > and (a,t) i— > h(a,t) be a C 2 function defined on 
(0,oo] x [a, b], which satisfies the heat equation 

. / / * / , m dh . d 2 h. 

V(<7,t)G(0,(7o] XK _( <7) t) =0 ._( (Tjt ). 

We assume that: 

(a) there are no t G [a, b] such that h(ao,t) = and ^(cr ,t) = 0, 

(b) there are no (a, t) G (0, ao] x [a, b] such that h(a, t) = and Vh(a, t) = 0. 

Then we have the following properties for the zero- crossings of h: 

(i) Global curves: If to G (a,b) is such that h(o~o,to) = 0, there exists 
o~q < o"o an d a maximal continuous path a i— > r to (cr) defined on (<7q~,<7o] such 
that r io (<To) = to and for all a £ (<Tq , <to] we have h(a, T to (<r)) = 0. Moreover, 
ifT to (a) stays within some compact set o/R for all a, then <Tq~ =0. 

(ii) Nonintersecting curves: If to ^ to is another point in (a,b) such that 
h(ao, to) = 0, then for all a G (0, ao] we have F to (a) ^ T t ~ (a) . 

(iii) Local description of the curves: If (a±,ti) G (0,a"o] x R is such that 
h(a\,ti) = then there exist a C 1 function rj defined on a neighborhood of 
a\ and such that h(a,n(a)) = in this neighborhood of a\, or a C 1 func- 
tion £ defined on a neighborhood of t\ and such that h(£(t),t) = in this 
neighborhood of t\, and moreover, if £'(t\) = 0, then £"(ti) < (it is a local 
maximum). 

The properties stated in Lemma 3 are illustrated on Figure 2, where the 
different types of curves formed by the set of points {{t,a) G R 2 ;/t(<r, t) = 0} 
are shown for some h satisfying the heat equation. 
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Fig. 2. Curves ofh(a,t) = for some h satisfying the heat equation, in the (t,a) domain; 
here t is along the horizontal axis and a is along the vertical one. According to Lemma 3, 
the zeros-crossings of h are a set of nonintersecting curves, that are locally else functions 
of a or functions of t with no local minima. 

Let us consider again the shot noise process X XjfT . We now give the main 
result for the number of local extrema of X\ a as a function of a. The 
intensity A is assumed to be fixed. 

Theorem 4. Let <r > and a<b. Then, 

P(3(cr, t) G (0, cr ) x [a, b] such that X Xa (t) = and VX Xa (t) = 0) = 0. 

Moreover, if we assume that for all < o\ < o"o 

E(#{cr G [cri, cr ] such that X' Xa (0) = 0}) < +oo, 

then the function a i— > p(X, a), which gives the mean number of local extrema 
of X\^ per unit length, is decreasing and it has the limit 2A as a goes to 0. 

Proof. Let us denote Y(a,t) :=X' Xa (t) for all (a,t) G (0,+oo) x R. 
We first check that the assumptions (a) and (b) of Lemma 3 are satisfied 
almost surely for Y. Assumption (a) is already given by Proposition 9. For 
assumption (b), we first notice that since Y(a,t) satisfies the heat equation, 
we have 

{Y(a, t) = and VY(a, t) = 0} 

= {Y{a,t) = and Y'(a,t) = and Y"(a,t) = 0}. 

Then a slight modification of the proof of Proposition 1, using the second- 
order Taylor formula in (10), allows us to conclude that P(3(<r, t) G (0,ao) x 
[a, b] such that Y (a, t) = and VY(<7, t) = 0) = 0, using the same integrabil- 
ity bound for the characteristic function of Y(a, t) as the one obtained in the 
proof of Proposition 9 [and considering first (a,t) G (cri,cro) f° r °~\ > 0, and 
conditioning by {77 1 > 3}]. This also proves the first part of the theorem. 
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Let < o"i < do be fixed. By assumption, we have E(#{cr G [01,00] such 
that X' Xa (0) = 0}) < +oo. Notice that by stationarity this expected value is 
independent of the value of t (taken as above). Let T > and let us consider 
the zeros of Y(a,t) = X' Xa (t) for (a,t) G [01,00] x [0,T]. Let t G [0,T] be 
such that Y(ao,to) = 0. By Lemma 3, there is a continuous path a i— > Tt (a) 
that will "cross the left or right boundary of the domain," that is, be such 
that there exists a G such that rt (0) = or T, or will be defined until a\ 
and such that T to (ai) G [0, T]. We thus have 

p{vo, [0,T]) < 2E(#{a G [01,00] such that X'^(0) = 0}) [0,T]). 

Dividing both sides by T and letting T go to infinity then shows that p{(Tq) < 
p(o"i). Thus the function a \- > p(X,a) is decreasing. 

To find the limit of p(X,a) as a goes to (that exists thanks to the 
bound of Proposition 9), instead of looking at the local extrema of X\ t(T in 
[0, 1], we will only look at the local maxima (which are the down-crossings 
of by the derivative) in [0, 11. Let D x > (0, [0, 11) be the random variable 

that counts these local maxima, and let p~{X,a) = E,(D X ' (0, [0, 1])). By 

stationarity of X\ t(T (t) and because between any two local maxima, there 
is a local minima, we have that p~(X,a) = ^p(X,a). Now, we introduce 
"barriers" in the following way: let E ao be the event "there are no points of 
the Poisson point process in the intervals [— 2a"o, 2o"o] and [1 — 2<7o, 1 + 2o"o]." 
If we assume that E ao holds, then X Xa {t) > for all t in [— 00,00] U [1 — 
00, 1 + 00] and all < 00, and therefore there are no local maxima of X\^ 
in these intervals. Then by Lemma 3, we can follow all the local maxima 
of X\ a in [0, 1] from = 00 down to = 0. Thus 1— > D x > (0, [0, 1])1e ct . 

A,er 

is a decreasing function of a for a < gq. Moreover, we can also check that 
the set of local maxima of X\ >(7 (t) in [0,1] converges, as goes to 0, to 
the set of points of the Poisson process in [0, 1]. This implies, in particular, 
that D x > (0, [0, 1]) goes to #{rj G [0, 1]} as goes to 0. Thus by monotone 

A, o" 

convergence, it implies that p~(X, a\E ao ) goes to E(#{rj G [0, lJljE 1 ^). Since 
the sequence of events E ao is an increasing sequence of events as 00 decreases 
to 0, we finally get 

hmp-(A,0)= hm E(#{r 4 G[0,l]}|K o )=E(#{^G[0,l]})=A. n 

Thus, under the assumption that E(#{<7 G [01,00] such that ^ CT (0) = 
0}) < +00 for all < 01 < 00, Theorem 4 asserts that the function a 1— > 
p(X, 0) is a decreasing function with limit 2A when a — > 0. This fact was em- 
pirically observed on Figure 1, and is also illustrated on Figure 3 where we 
"follow" the local extrema as evolves. Now, these properties can be trans- 
lated, using the scaling relations of Lemma 2, into the following properties 
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Three processes t -> X (t) for respectively a = 0.1 ; 0.3 and 0.8 




Fig. 3. Top: three processes t 1-> X\, a (t) obtained from the same Poisson point process of 
intensity A = 2 and for a Gaussian kernel of respective width a = 0.1; 0.3 and 0.8. Bottom: 
evolution of the local extrema of f H X\ :Cr (t) as o goes from to 1. The three values 
a — 0.1; 0.3 and 0.8 are plotted as dotted line. They indicate the local extrema of the three 
processes above. 



on A i—)- p(X, a): 

Vc > 1, p(cX, a) < cp(X, cr); 



p(A, cr) < 2A and 



p(X,cr) 



1. 



2A a^o 

This shows the second asymptotic linear regime observed for small values of 
the intensity A. 

APPENDIX 

A.l. Stationary phase estimate for oscillatory integrals. 

Proposition 10 (Stationary phase estimate for oscillatory integrals). 
Let a <b and let ip be a function of class C 2 defined on [a, b] . Assume that ip' 
and ip" cannot simultaneously vanish on [a, b] and denote m = 
min sg [ a b ] yj tp'{s) 2 + Lp"{s) 2 > 0. Let us also assume that uq = #{s € [a, b] s.t. 
tp"(s) = 0} < +oo. Then 



Vu G E s.t. \u\ > 



rn 



f 

■J a 



< 



8\/2(2no + 1) 



m\u\ 



Now, let (fi and if2 be two functions of class C 3 defined on [a,b]. Assume 
that the derivatives of these functions are linearly independent, in the sense 
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that for all s G [a, 61, the matrix &(s) = i^yffl ^f^fl) is invertible. Denote 
J w Vi(s) v 2 ( s ) 

sG[a,fe] II^K 5 ) -1 !!" 1 > 0> where || • || is the matricial norm induced by 

the Euclidean one. Assume, moreover, that there exists uq < +oo such that 

#{s£[q,6] s.t. det($'(s)) = 0} < n , w/iere 



??? 



mm 



8V^(2n + 3) 



. T/ien 



V(ii,u)eIR 2 s.t \f- 



u A + v z > — 
m 



Jwpi(s)+ivip 2 (s) ^ g 



< 



V m\fu l + w 2 



Proof. For the first part of the proposition, by assumption, [a, 6] is the 
union of the three compact sets 

{s £ [a,b];\tp"\ > m/2}, {s G [a, b]; \ip'\ > m/2 and ip" > 0} and 
{s G [a, 6]; \ip'\ > m/2 and y?" < 0}. 

Therefore there exists 1 < n < 2uq + 1 and a subdivision (aj)o<i< n of [a, b] 
such that [aj_i, aj] is included in one of the previous subsets for any 1 <i<n. 



page 332, 



otherwise, 



] C {s G [o,6];| 




s)\ > m/2}, according t 








/ e iuv(s) ds 




I""* e iu(m/2)(2<p(s)/m) dg 


J 04-1 







< 



v/2 



^m\u\ 



ra,i 

/ 



< 



m\u\ 



The result follows from summing up these n integrals. 

For the second part of the proposition, we use polar coordinates, and write 
(u,v) = (rcos#,rsin#). For 9 G [0,2-/r), let fg be the function defined on 

[a,b] by tp 9 {s) = Ms) cosO + ip 2 (s) sm6. Then = and thus 

1 = ||*(a) -1 (^)||. This implies that for all s G [a,b], vV (s) 2 + </ '(s) 2 > 

l/||$(s) _1 || > m. Moreover, thanks to Rolle's theorem, the number of points 
s G [a, b] such that ip' e ' (s) = is bounded by one plus the number of s G [a, b] 
such that y>'{(s)ip2 \s) — <Pi{s)(p 2 ' (s) = 0, that is, by 1 + no- Thus, we can 
apply the result of the first part of the proposition to each function ipg and 
the obtained bound will depend only on m, no and r = \fu 2 + v 2 . □ 



A. 2. Proof of Proposition 8. For k > and I > integers, let us denote 
m ki = f \9(s)\ h \g'{s)\ l ds. We will also simply denote mo = m2o = / g{s) 2 ds 
and m-2 = mo2 = / g'(s) 2 ds. 
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Let ip\(u, v) denote the joint characteristic function of (Z\(t), Z' x (t)), then 
tp x (u,v)= E( e ! ( u /^+ i W^) Jf i)e-™ v/X / a 

= exp^A J f e <u/VX)9(s)+i(v/VX) 9 '( s ) _ J _ i-j=g(s)) ds) . 

We now use the fact j g' = 0, and we thus have ip x (u,v) = ex.p(H x (u,v)) 
where 

H x (u,v) = \ J (e i ^^(s)+i(v/^)g\s)_ 1 _ i ^ g ^_ i v,^\ ^ 

We need to notice that 

V(u,v) G R 2 \Mu,v)\ = \exp(H x (u,v))\ = \E{e luZ * +ivZ '*)\ < 1. 
In the following, we will also need these simple bounds: 



(28) 



e ix - 1 - ix + — 



I X I 

< -gj- and 



VzeC |e z - 1| < |z|e |z| . 
We first estimate H x (u,0). We have 

H x (u,0) = \ J hw^Ms)_ 1 _ i JL g{s) \ ds = _}_ u 2 mo + Kx{u)j 

where K x {u) = \J( e i( - u /^ 9 ^ - 1 - i^g(s) + \^g 2 {s))ds. Then, thanks 
to the simple bounds (28), we get 

|-^a( u )| < ■ — • ? ^- 3 ° and consequently 
6v A 

\ e H x {u,Q) _ e -(l/2)« 2 m | < H!^30 e -(l/2)u 2 mo e |u| 3 m3o/(6v / A) 

~ 6^ 

We then estimate H x (u,v) — H\(u,0), 

H x (u,v) - H x (u,0) = A / ^/VX)g{s)+i{v/^)g'{s) _ e i(u/V\)g(s)^ ds 



= -y / g'(s) 2 e^/^ s Us + F x (u,v), 

where F\ (u, v) = A / e i(«/VA)tf(«) ( e i(«/vW») _ i _ i-^</(s) + g 5 '( s ) 2 )ds. 
And again, thanks to the simple bounds (28), we get |-F\(u,t;)| < ^J^ 3 . 
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This implies that 

\ e H x (u,v)-Hx(u,0) _ e -f 2 / 2 / g'{s) 2 e i( - u /^3M ds , 

< | e -f 2 /2 / a 1 , ( s )2 e i C«/vO:) 9 ( s ) ds ^ ,f x (u,v) _ ji 

< b| 3 m 3 c _ t , 2/2 r g '( s )2 C o S ((«/VA)g(g))rfs4>| 3 mo 3 /(6^) 

~ 6^A 

Let us now compute Cz x (u, [0, 1]). By Proposition 1, we know that 

-ttC Zx {u,[0,1})= / -z(i/>x(u,v) + il>x('u,-v)-2il)x{u,0))dv. 
Jo v 

Let y > be a real number. We split the integral above in two parts, and 
write it as the sum of the integral between and V, and of the integral 
between V and +oo. Since for all (u,v), we have \i/jx(u,v)\ < 1, we get 



+oo 2 

-o(tpx(u,v) +Tp x (u,-v) - 2ipx(u,0))dv 

V v 



■+oo j 4 

<4/ ^ dv = — . 



On the other hand, let Iv(u) denote the integral between and V. We have 

Iv (u) = f V \ e H X (u,0) {e H x (u,v)-H x (u,0) + e H x (u,-v)-H x (u,0) _ y ^ 



JO u 

We then decompose this into 

Iy(u)= / _ e H X {u,0), e H x {u,v)-H x {u,0) + e H x {u-v)-H x {u,0) 

Jo v 2 

+ f ]_ e H X (u,0) ( 2e - 1 ' 2 /2/ 9 '(s) 2 e^/v^) 9 (s) ^ _ 2e _(„2 /2)m2 ^ 
rV i 

+ / _L( e ^A(«,0) _ e -(l/2)« 2 m + e -(l/2)« 2 mo\(2 e -(^/2)r7i2 _ 2 ) ^ 

Using the bounds we computed above, we get that 

rV „-<v 2 /2)m2 1 

I v (u)-2e-^> 2 ^ e \ ~ l dv 
. n v z 



< 2 



V vm 03 ^_ v 2 /2 f ql[s? cos((u /^x )q(s)) ds +\v\ 3 m 3/((>J\) dy 
6\f\ ' 



rV 

+ 2 



l| e -i, 2 /2/ 9 '(s)V("/^) 9 (s) ds _ 2e -(v 2 /2)m 2 \ dv 







,fV i „-(« 2 /2)m 2 

+ 2 | e H A Ko)_ e -(i/2)« 2 m o| / 1 - e V ; dv . 



V 2 
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Let Jy\u), for n = 1,2,3, respectively, denote the three terms above. To 
give an upper bound for Jy (it), we will need the following basic inequality: 



vm 03 c -v 2 m2/2+(v 2 /2)(u 2 m22/(2X))+\v\ 3 m 03 /(6VX) dv 



Vx G R, cos(x) > 1 - \. This gives us the bound 

4\u)<2 [ 
Jo 

For the second term, we use 

| e -« 2 /2/g'(s)V( u /VA)s(s)(2 s _ 2e -(« 2 /2)m 2 | 

<• e -(t ; 2 /2)m2| e -(i) 2 /2)/5'(s) 2 (e i ( 11 ' / ^)9( a )-l)ds _ 



< e -(«72)m 2 



ff '( s )2( e i(«/VA) fl ( S )_ 1 ) ds 



J(f 2 /2) / 9 '(s) 2 (e ! ( 1 '/^)9( s )-l)ds| 



But \jg l {s) 2 {e i ^l^9^-l)ds\<jg\s) 2 ^g{s)ds = ^m l2 and thus 

4 2 )( U ) < H mi2 [ V e -(v 2 /2 )m2 + (vV2) { \u\/VX) mi2dv _ 

v A Jo 

For the third term, we use an integration by parts to obtain that 

! _ e -(, 2 /2)m 2 e -(V 2 /2 )m2 _ ! ,V : 



/ ^— ^ dv = - -+/ m 2 e-^ 2 /2)m 2 ^< 

Jo v * Jo 

which gives 

4\u) < V2^M!^ e -a/2)u 2 m +| U |3m 3 o/(6^)_ 



vrm 2 , 



6^A 

Moreover, we also have 

V 1 _ e -(« 2 /2)m 2 



/•v 1 _ e -^-/2;m 2 1_ e -(V 2 /2)m 2 /-+oo 

/ _ cfo - v/2i^ < - — - — — + / m 2 e-( v dv 

Jo v V Jv 

2 



< 
~ V 

The partial conclusion of all these estimates is that 
| irC^ x (u, [0, 1] ) - V2^2~e~ mou2 / 2 \ 

<y + ^ + 4\u) + J { y\u) + 4\u). 

We now have to choose V in an appropriate way. The choice of V will be 
given by the bound on Jy (it). Assume in the following that it satisfies the 
condition (Z71) given by " 2 ^ 22 < and let us set 

_ 3\/Am 2 
4m 3 
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Then for all v G [0, V], + + < and thus 

4\u)<^ [ V ve-^dv< 2m ° 3 / _ . 

3v A jo 3m2vA 

For the term J y (u), we notice that if u satisfies the condition (U2) given 
by 77=mi2 < ^r, then for all V > 0, we can bound Jy (u) by 

J«(«) < ^=m 12 f V e -^l^ dv < M 

VA Jo VA V m 2 

Finally, for the third term, we have that if u satisfies the condition (£73) given 

u\m S o ^ l ,i 7 (3) 

3v /y S 2 m o, then J y 



by 1"^° < imo, then Jy^(it) can be bounded, independently of 1/, by 



t(3)/ \^ A5 1«| 3 ^30 -(l/4)« 2 mo . /o 2|«|m 30 _i 

6vA 3movA 

because of the fact that for all x > 0, then xe~ x < e" 1 . 

The final conclusion of all these computations is that if we set a\ 

min (7S' & then for a11 ^ and A > we have 

|u|<ai>/A => |^ A (n,[0,l])- V2i^2"e- mou2 / 2 | <^= + ^, 

v A v A 

where a 2 = ^±^a and a 3 = m 12 ^ + 



3m 2 ^""^y m 2 1 3m 

A. 3. Proof of Lemma 3. The proof of this lemma relies upon the implicit 
function theorem. Let us start with the proof of (i); let (o"0)*o) be a point 
such that h(ao,to) = 0. By Assumption (a), we have that ^r(co,to) 7^ 0. 
Then, thanks to the implicit function theorem, there exist two open inter- 
vals I = (ctq~,o"q) and J = (ip >*o") containing, respectively, ctq and to, and 
a C 1 function rj:I — > J such that T]{cfq) = to and V(<7, i) G I X J, h(a,t) = 
^ t = rj(a). Let us now denote 7? = r to . We need to prove that we can 
take o"q~ = when Tt remains bounded. Assume we cannot; the maximal 
interval on which Tt is defined is (<7^,<7q~) with > 0. By assumption, 
there is an Mo > such that for all a G (<Tq~,<7q~), then |r to (o")| < Mo. We 
can thus find a subsequence (cr^) converging to Oq as fc goes to infinity and 
a point t\ G [—Mo, Mo] such that r to (<7fc) goes to ti as fe goes to infinity. 
By continuity of h, we have h{a^ ,t\) = 0. Now, we also have ^(<7q~, ii) = 0. 
Indeed, if it were 7^ 0, we could again apply the implicit function theorem in 
the same way at the point {a^ ,t\), and get a contradiction with the max- 
imality of I = {0^,0$ ). Then, by Assumption (b), we have ^{o^ ,t\) / 0. 
We can again apply the implicit function theorem, and we thus obtain 
that there exist two open intervals I\ = (a^,af) and J\ = (t^,tf ) con- 
taining, respectively, and t±, and a C 1 function £: J\ — > I\ such that 
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£(ii) = <T and V(cx, t) 6 I\ x J\, h(a,t) = 44> cr = £(t). Moreover, we can 
compute the derivatives of £ at ti. We start from the implicit definition of 
£: h(£(t),t) = 0. By differentiation, we get £'(i)§(£(i),t) + §(£(i),t) = 0. 
Taking the value at t = t±, we get £'(ii) = 0. We can again differentiate, and 
find £"(t)^(£(i),t) + £'(^0(^( t ) )t ) + 2£'(i)S(£(t),i) + 0(£(t),t)=O. 
Taking again the value at t = ti, we get 

Thus it shows that £ has a strict local maximum at t%; there exist a neigh- 
borhood U\ of <7q~ = £(ti) and a neighborhood Vi of ti such that for all 
points in U\ x Vi, then h(a,t) = implies cr = £(t) < £(ti) = ctq~, which is 
in contradiction with the definition of Ft on (<7q~, Cg - )- This ends the proof 
of (i), and also of (hi). 

For (ii), assume that to an d to are two points such that ft,(o"o,fo) = 
h(ao, to) = and such that there exists o\ < a$ such that T to (<Ti) = T t ~ (<7i) = 

t\. Then, if ^r(ci, ti) 7^ 0, the implicit function theorem implies that Tt (cr) = 
r t ~ (cr) for all a G [a±, o"o] and in particular to = But now, if ^r(ci, ti) = 0, 

then, as above, this implies that |^(cji,ti) 7^ and using again the implicit 
function theorem, this would be in contradiction with the fact Tt (cr) is 
defined for a G [a\ , Co] . 
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